CENTRAL LIMIT THEOREMS FOR THE LARGE-SPIN 
ASYMPTOTICS OF QUANTUM SPINS 



TOM MICHOEL AND BRUNO NACHTERGAELE 

Abstract. We use a generalized form of Dyson's spin wave formalism to prove several 
central limit theorems for the large-spin asymptotics of quantum spins in a coherent 
state. 



1. Introduction 

In statistical mechanics, the thermodynamic limit of an infinite number of interacting 
particles in the continuum or on a lattice can be taken rigorously using the laws of 
probability. In a first approximation we are usually interested in the behavior of intensive 
observables, i.e., observables that grow proportionally to the total number of particles. 
Taking their thermodynamic limit corresponds to the Law of Large Numbers (LLN) in 
probability theory. Introducing quantum mechanics at the microscopic level does not 
change much at the macroscopic level: intensive observables still behave classically in 
the thermodynamic limit, or, the LLN for quantum systems, in particular, involves a 
description of the system in terms of classical, i.e., commuting, variables. 

The next logical step is to study fluctuations of intensive observables around their mean. 
The corresponding law in probability theory is the Central Limit Theorem (CLT). Here, 
quantum effects can survive in the limit of an infinite number of particles, and fluctuations 
often behave non-classically and have to be modeled using non-commuting variables. 
Typically, this happens in the presence of a spontaneously broken continuous symmetry, 
and the macroscopic quantum fluctuations present in such a system are the well-known 
Goldstone bosons pQ. A general theory of Goldstone bosons using non-commutative 
central limit theorems was presented in [TT] . 

The study of non-commutative central limit theorems becomes essential in this context, 
and a general theory was developed in [HIGH] ■ The main result is that the macroscopic 
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quantum fluctuations can be identified with a representation of the Canonical Commuta- 
tion Relations (CCR) in a quasi-free state which is determined by the correlations in the 
microscopic state. The non- commutative central limit theorem can be considered as the 
first quantum correction to the LLN for these systems. 

In addition to the thermodynamic limit (i.e., the large- N limit), also the so-called clas- 
sical limit results in an asymtotic description of the system by classical random variables. 
A well-known example is the classical limit of quantum spins, which was treated rigor- 
ously in great detail by Lieb ^5]- This is the limit of infinite spin: the quantum spin 
operators are normalized such that, in a suitable sense, they converge to classical spin 
variables with values in the unit sphere in R 3 . 

As an example of such a result we mention the classical limit of the free energy, obtained 
for a large class of models by Lieb In this case, the classical limit result tells us 

that after rescaling each spin operator in the Hamiltonian, the partition function of the 
quantum model will converge to the partition function of the corresponding classical 
model. One may also consider the limit of large spin in a sequence of states, and compute 
expectations, in which case one is interested in the classical limit as a result about the 
distribution of the rescaled spin operators considered as random variables. In this present 
work we are interested in the latter situation. 

In Lieb's treatment of the classical limit an important role is played by the so-called 
coherent states, which were studied in great detail in Among all quantum states, 
the coherent states are those that optimally approximate the idea of a spin pointing in 
a certain direction in space. Mathematically this is made precise by showing that these 
states carry no dispersion for the classically rescaled spin operators (see Proposition 12. 
Again, we see that the classical limit for spins, which is similar to a LLN in probability, 
results in classical random variables. The role of N is played by the magnitude of the 
spin, which we will denote by J. It is now reasonable to expect that in such states in 
which the 'intensive' (scaled) observables become dispersionless in the limit, that a central 
limit theorem should hold for the fluctuations of these observables around their expected 
mean. The subject of this paper is then the formulation and proof of various central limit 
theorems for these fluctuations, which are represented by operators on a Hilbert space. 
The class of states we consider are products of so-called coherent states for the spins. 

In Section El which also contains the mathematical setup including the definition and 
basic properties of the coherent states, we formulate a first non-commutative Central 
Limit Theorem (CLT) for large J using the techniques of Goderis, Verbeure, and Vets 
[8|ll(J]. In a rather different setup, a similar result was previously given by Michoel and 
Verbeure in [16] . The proof of this CLT and some auxiliary results are given in Appendix 

m 

The main purpose of this paper is to strengthen the results of Section EJ which could 
be called standard, by taking advantage of the additional structure that is present in 
models with a high dimensional representation of the SU(2) commutation relations. More 
precisely, we want to use Dyson's spin wave formalism (ii. which can be briefly described 
as follows. The 'all spin up' state is used as a reference state in spin Hilbert space. It plays 
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the role of the vacuum state in the sense that there is a formal analogy between lowering 
the spin in this state, and creating a boson particle in a Fock vacuum state. Dyson made 
this analogy precise by defining a unitary equivalence between the spin Hilbert space and 
a subspace of Fock space (the subspace with no more particles present than the size of 
the spin). Under this equivalence, the spin lowering operator, divided by the square root 
of the size of the spin, becomes the boson creation operator up to a correction that is a 
function of the number operator that formally goes to one in the large-spin limit. After 
the trivial observation that every coherent state is the 'all spin up' state for the spin in 
the defining direction of the coherent state, we see that Dyson's formalism can be used to 
study the coherent states of our situation. The spin wave formalism is discussed in detail 
in Section El 

The scaling with the square root of the spin in Dyson's equivalence between the spin 
lowering and the boson creation operator is precisely the scaling we need in the central 
limit theorem for the fluctuations of the spin operators. Hence, for all values of the spin, 
we can write the fluctuation operators as well-defined operators on one and the same 
Fock space, and it follows (see Section Hj), that the limit of infinite spin can be taken as a 
genuine operator limit. This is where the central limit theorem for spin fluctuations differs 
fundamentally from the known non-commutative central limit theorems and it makes it 
possible to prove stronger results. For instance, the usual central limit theorems prove 
convergence of the characteristic function of linear combinations of fluctuation operators. 
Using Dyson's spin wave formalism, we obtain convergence of the characteristic function 
of arbitrary polynomials of fluctuation operators. 

In the last section, Section we discuss some applications of this stronger version of 
the central limit theorem. 

The first application (in Section 15.1)) is to obtain a rigorous version of "bosonization" 
for quantum spin systems. This refers to the well-known technique in physics that the 
low-energy excitations, and hence the low-temperature behavior, of certain quantum spin 
systems can be well approximated using a boson or spin wave approximation. This idea 
can be made mathematically precise for those quantum spin systems which possess a 
coherent ground state. Indeed, using our central limit theorem we immediately obtain 
strong convergence of the spin Hamiltonian and the spin dynamics to a quasi-free boson 
Hamiltonian and dynamics. In [THj, we use this convergence to obtain the large-spin 
asymptotics of the energy spectrum of the anisotropic ferromagnetic Heisenberg chain, 
thus improving upon earlier results We also discuss another application of this 

convergence result, namely to obtain the large-spin asymptotics for the time evolution of 
the ground state when the dynamics is perturbed by a spin fluctuation. 

In a second application (Section 15.2)1 . we apply the theorem to the study of non- 
commuting fluctuation operators for N independent copies of quantum random variables. 
A central limit theorem for the characteristic function of polynomials of such fluctuation 
operators, analogous to the one we obtain in Section has appeared in the literature in 
the form of a conjecture by Kuperberg [T5] . The special case of this conjecture for expec- 
tations computed in a tracial state, was proved in (Hj and establishes an interesting result 
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about the distribution of the shape of a random word. This is an example of a quantum 
CLT used to prove results in classical probability theory, in this case generalizing earlier 
work of Johansson [11]. 

The situation with a tracial state can be regarded as intermediate between classical 
and quantum probability theory, because the cyclicity of the trace implies that although 
the microscopic variables do not commute, the limiting fluctuations are classical Gaussian 
random variables. We can apply our results to obtain the first example of a non-tracial 
state for which Kuperberg's conjecture holds. We show that a system of N independent 
spin- 1 particles, each in a coherent state, can be identified with one spin-^ particle in the 
corresponding coherent state. Therefore, our general results can be applied directly and 
we obtain a fully non-commutative system in which the conjecture is valid. 

2. Mathematical setup and preliminary results 

We consider quantum spin systems on a finite or infinite lattice £. At each site x G £ 
we have a spin-J degree of freedom (J G |N ), i.e., a (2 J + l)-dimensional irreducible 
representation of SU(2), and we denote with 5* the corresponding spin-J matrices, 

[S l x , S J y ] = i5 x , y £ijkSx 

S x -S x = (Sl) 2 + (S 2 x r + (Slf = J(J + l) 

We will also use the spin raising and lowering operators: and S~ , = ± iS%. 
The local Hilbert spaces associated to each finite subset A of the lattice are therefore 

and the algebra of observables is 

21^ = 0(^12^+1(0))^ 

xeA 

For A' C A, 21j ; a' can be considered as a subalgebra of SI^a in a natural way. With this 
in mind, the union 

2tj,ioc = |J 2l A ,j 

is the algebra of local observables, and its closure is the quasi-local algebra of observables 
2lj of the spin system on £. For a finite lattice £ this is of course just the local algebra 
specified before. 

As a basis for the local Hilbert spaces we take the standard basis which is diagonal for 
the S x operators: 

Sl\m x ) = m x \m x ), m x = -J, -J + 1, . . . , J 
Sx\ m x) = a/J(J+ 1) - m x (m x ± l)\m x ± 1) 
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A coherent spin state j2JEi] at site x is specified by a unit vector in R 3 (i.e., a classical 
spin) u x = (9 X , (p x ), < 9 X < 7T, < (p x < 2tt, and defined by 

\(9 x ^ x )) = e^^- s ^-^\J) 

J / 9/ \V 2 

= V L (cosi^) J+m ^(sini^) J -^e i(J - m ^^ \m x ) 

The states of the spin system that we will consider are tensor products of coherent 
states: 

u{A) = <8> ((6 x ,(p x )\A\{6 x ,(p x )) 

which is well-defined for all A e 21 j. The collection of unit vectors {u x = (9 x ,(p x )} xe £ 
defining the state is completely arbitrary. The GNS Hilbert space of u is the incomplete 
tensor product Hilbert space 



Hj= |J ( [ ® C2J+1 ] ® [ ® fi 



Ac£ 



where f2 y is short-hand for Q y = \ (9 y , <p y )), and the generating vector is 

Q j = ® Q x 

Because of the simplicity of the GNS representation we will not need to distinguish between 
an observable and its representative in B(Hj). 

A particular type of quasi-local observables that we are interested in are those of the 
form 

^2v x -S x , v x eR 3 (1) 

In order that functions, such as polynomials and trigonometric functions, of these ob- 
servables are well-defined, we need summability conditions on the set of v x . Consider 
v = {v x G IR d }a; G £, and define the p-norm (p > 1) by 

where \v x \ is the length of v x . The space of all Revalued summable sequences (i.e., \v\i < 
oo) on £ will be denoted ^(£), and the space of all Revalued square summable sequences 
will be denoted ^(£). For d = 2 we get the standard complex valued (square) summable 
sequences, and for these spaces the subscript d is omitted. Note that if \v\ Po < oo for 
some po > 1, then for p > p , \v\ p < \v\ Po < oo as well. 

In a number of generating functions are derived for coherent states. The most 

important one for us is the characteristic function for observables of the type (|T|): 

w ( e *E.n.-S.) =n{cos(^,|)+^^^sin(i|^|)} 2J (2) 



\ v \p ■= 2^1^ 
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It follows immediately that 



uj 



{jjVx ■ S x ) = J^Vx 



u. 



The classical limit of a quantum spin system is obtained by normalizing each spin 
operator by J, and then taking the limit J to oo. The quantum spin variables then 
converge to classical spins taking values in the unit sphere in R 3 . There are various 
precise mathematical statements that can express this. The most common are results 
formulated as convergence of the free energy of a quantum spin system to the free energy 
of the corresponding classical spin system In the present paper, we don not specify a 
Hamiltonian, but consider the characteristic functions of the spin variables in states that 
are products of coherent states. The following result, proved in Appendix |X] as an easy 
consequence of eq. (J2J), is a law of large numbers for the quantum spins. 

Proposition 2.1 (Classical limit). For v £ ^>(£) and J £ |Nq, we have 



■( 



uj [e-J 



7 J2 X v ^-Si 



) 



iJ2 x Vx-Ux 



< 7 exp 



|u|i + 2exp(|u|i) 



More general results for products of characterisitc functions, analogous to Theorem 12.31 
below, can also be derived for the classical limit, but the above proposition is sufficient 
for our needs. 

The next step is to study the fluctuations of these observables around their mean and 
we define fluctuation observables by 



Fj{v) = J j^[vx ■ S x - u(v x ■ S a 
xez 



(3) 



Our first result here is about the characteristic function of the fluctuation observables in 
a coherent state: 

Proposition 2.2. For v £ and J £ |N , we have 



where 



e a 1,2 



< 



J 1 /' 



: b(v) 



b(v) 



exp 



2 1/ >| 2 + 2 1/2 exp(2 1/ >| 2 ) 



(4) 



The vector v £ 

is defined as follows. At each site x we start with the same standard basis fel,e^,e 3 ) of 



! (£) in the proposition whose length determines the variance of Fj(v) 

■'xi^xi' 



R . Through the state uj, we are given a unit vector u x 



'XI 



ip x ) where the spherical 



coordinates are given w.r.t. the standard basis. We can define a new basis of R 3 (now site 
dependent) by rotating e 3 to u x . The rotated e\ and e\ then span the plane orthogonal 
to u x which can be identified with the tangent plane to the unit sphere at u x . We define 
v x £ R 2 to be the projection of v x onto this plane. 
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To see that this is really the variance, we can compute directly from eq. (J2J that 

lo(Fj(v)Fj(w)) = S ^[v x - w x - (v x ■ u x )(w x ■ u x ) + i(v x x w x ) ■ u x ] (5) 

and it follows that indeed 

uj{Fj{v) 2 ) = \v\l 
The proof of Proposition 12.21 is given in Appendix 1X1 as well. 

We can continue along the lines of [EJEI1 and show that in the limit J —>■ oo the system 
of fluctuations is given by a representation of the canonical commutation relations (CCR). 
The most general theorem in this context is: 



Theorem 2.3. For n G N , fi, . . . , v n G and J G |N , 

Tl 71 /* 71 71 — 1 71 n 

Kik^O -*(n»m>)| * +e«("*- e «*)} <«) 



where it is understood that the second sum is zero for n = 1, b(v) is given in (@J) and 
a(v, w) is given by 

a(v, w) = -\v \2\w\2Hv | 2 + \w\ 2 ) + v^exp ^|v| 2 |w| 2 + exp(|v| 2 |w| 2 ) 

In this theorem, the Vj G £ 2 (£) are again defined by projecting onto the tangent 
planes at the different u x \ the W(v) are the Weyl operators generating the CCR algebra 
CCR(£ 2 (£), a), and a is the symplectic form associated with the standard inner product 
in £ 2 (£), i.e., 

a(v, w) = 23xn (v,w) = 23m S y^ j v x w x 

and Cj is the quasi-free Fock state defined by 

u(W(v)) = e-^'V = e-sWi 
The Weyl operators satisfy the well-known canonical commutation relations: 

W(v)W(w) = e'^^Wiv + w) 

The proof of Theorem I2.3I requires a Baker-Campbell-Hausdorff-type formula (BCH 
formula) that shows that the e lFj ^ approximate these commutation relations in a suitable 
sense, and then proceeds from Proposition I2.2I through a standard induction argument. 
For convenience of the reader, complete proofs are included in Appendix IX! 

The proofs of both Proposition I2. II and I2.2I and Theorem 12 . 31 relv directly on the use of 
the generating function (J2J). In order to prove convergence of moments or convergence of 
the characteristic function of polynomials of the Fj(v), Dyson's spin wave formalism will 
prove much more convenient. 
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3. Dyson's spin wave formalism for coherent states 

In his famous papers 0111], Dyson introduces a formalism for studying rigorously 
bosonization in the Heisenberg ferromagnet. The main idea is to take the ferromag- 
netic ground state Cg^es | J) as a reference state, and to identify the lowering of spins in 
this state with the creation of boson particles in a Fock state. This allows the Hamiltonian 
to be written as an operator on Fock space, and it is argued that the leading terms at low 
temperature are the ones linear and quadratic in the creation and annihilation operators, 
thus obtaining an exactly solvable system. This approximation corresponds to taking a 
large J limit much as we want to do here. A nice exposition of Dyson's formalism in a 
more modern language is in [22], and a rigorous theorem about convergence of the free 
energy within this formalism is in [3]. 

Upon inspection it is clear that the same formalism can also be used if the reference 
state is a general product state of coherent states instead of the purely ferromagnetic state 
® X \J). This was used in [TH] to prove that the low-energy excitations of interface ground 
states of the 1-dimensional XXZ chain are given, in the large J limit, by a quadratic boson 
Hamiltonian describing particles hopping on the lattice under the influence of an external 
potential centred around the interface. In this paper we use the formalism in an analogous 
way to prove several central limit theorems for the fluctuation observables Fj(v). 

Recall the new basis we defined for every site x using the unit vectors u x that determine 
the state ui (see after Proposition 12.2(1 . More precisely, this rotated basis is given by 

f x = cos 9 X cos ip x e\ + cos 9 X sin (p x e x — sin 9 X e x 
fx = ~sm (p x el + cos(p x e 2 x 

f x = u x = sin 6 X cos <p x e x + sin 9 X sin ip x e x + cos 9 X e x 
The spin operators can be rotated likewise, and we find (using, e.g., eqs. (3.9) of [2]) 



S x 


= U X S X U X 


= fl 


' S x 


= cos 9 X cos ip x Si + cos 9 X sin Lp x S x — 


sin 9 X S Z X 


s~ 


= U X S X U* 


= fl 


' S x 


= - sin ip x Si + cos if x Si 




s x 


= U X S X U* 


= fl 


■ S x 


= u x = sin 9 X cos ip x S x + sin 9 X sin ip x 


S 2 x + cos9 x S 3 X 



where 

U = ^{S^e^-S+e-*^) 

The main observation is that in this new basis the state <g> x \ (9 X , (p x )) becomes the usual 
reference state ® x \ + J) for the spin wave formalism: 

S 3 x \(9 x ,<p x )) = SlU x \J) = U x Sl\J) = J\(9 x ,cp x )) 

The rotated spin raising and lowering operators are 



&x ~ iS x 
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or 

S+ = - sm0 x S 3 x + cos 2 {\d x )e- llp *St - sm 2 {\d x )e llp * S~ 
S+ = - sm9 x S 3 x - sm 2 {\O x )e ilp * S+ + cos 2 {\6 x )e- i{f> * S x 

According to our previous notation we denote the components of a vector v x £ R 3 w.r.t. 
the standard basis with v l x , and w.r.t. the rotated basis with v l x . We find 

1 x, 1 



v x -S x = -v x S+ + -v^S x + v 3 x S 3 x 



where v x = v x ± iv 2 , and 



*(») - E/."(J^ + (§r«S<3 - J) 



where we used u(v x • S x ) = Jv x ■ u x = Jv x . 
Let 



n = {n x £ N} xe£ 
Afj = jn | Vrr: n x < 2 J, ^ < oo| 



i 



where Q,j = ® x \ (9 X , (p x )) is the GNS vector for the state u. The set {(pa \ n £ A/j} is an 
orthonormal basis for Tij, which will prove to be very useful. 

Recall that Co is the quasi-free Fock state on CCR(£ 2 (£), a). Its GNS Hilbert space is 
the usual Fock space T with a vacuum vector Q, and creation and annihilation operators 
a* and a x : 

\Oj X) CLy\ 5 x y 

a x fl = 

The representative of W(v) is the usual Weyl operator 

W(v) = e iF ^ = e i ^" c ^ a * x+ ^ ax 
where, if v x = (v x -,v 2 ) £ R 2 , v x £ C are defined as above. The unbounded operator 

F(v) = Y.K<+^ a *} ( 7 ) 

is well-defined for v £ £ 2 (£) and is called the boson field operator. 
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Now let 

J\f — | n | n x < oo | 

X 



then we have an orthonormal basis {(p^ | n G A/"} of T . Identifying (p^ with cp^, it is clear 
that the GNS Hilbert spaces, Hj, J G |N , can be identified with a nested sequence of 
subspaces of J 7 , defined for each J as the linear span of all vectors (p^ with n x < 2J . 
More precisely, we use the projections P UjX on T which project onto the first 2n boson 
states at site x, i.e., on the states (p^ with < n x < 2n, and denote P n = Y[ x F n,x, and 
find 

Hj = PjT 

where = means unitarily equivalent. 

Under this equivalence, we find that the spin operators are given by (see |22] f° r more 
details) 

= Pja* x gj(xr\ j0— = gj(x)^a x Pj, J - S 3 X = Pja* x a x Pj 

where 

9j(x) = 9j(a* x a x ) 

and 



9 An) 

Hence 



I - jjn n<2J 
n>2J 



= Pj{j2^<9Axr 2 + v x gj(x)^a x - (^) 1/2 fel^ 



(8) 

If we let J — > oo, we have, for fixed n, gj(n) — > 1 and formally Fj(v) — > F(v), the boson 
field operator of eq. (J7J). To make this into a mathematically precise statement is the 
subject of the next section. 

4. Operator convergence of fluctuation operators 
For any n G N, let T> n denote the linear span of the vectors 

F(wi) . . . F(w n )n : w h . . . , w n G £ 2 (£) 

and define V by 
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i.e., T> is the linear span of the vectors ip^, n G N '. The set T> consists of entire analytic 
vectors for F(v), v G £ 2 (£) (Theorem 4.6 of [20.), i.e., ip G T such that ip is in the domain 
of F{y) k for every k G N and 

j.k 

E^II F (^ll <o ° ( t>0 ) 

fc>0 

and therefore T> is dense in T . It also follows that T> is a core for all F(v), v G £ 2 (£) as 
well as for the creation and annihilation operators a* (if) and a(v), 

a*(v) = [a(v)\* = J2v+a* x 

Clearly, any function of the number operators, n x = a x a x , such as Pj, leaves the spaces 
T> n invariant and hence: 

Lemma 4.1. For all n G INI, v±, . . . , v n G ^j(£), and J G |N suc/j £/ja£ 2 J > n, we have 

Fj{v 1 )...F J {v n )tteV n 
Recall the following lemma, Lemma 4.5 of |2"U] : 
Lemma 4.2. For all n G INI, tp n G D n and 5 G ^ 2 (£), we nave 

11^(13)^11 <2|i5| 2 (n + l) 1 / 3 ||Vn|| 
The analogous result for the Fj(v) is: 
Lemma 4.3. For all n G N, G f G ^(-C), arad J G |N smc/i i/iai 2J > n, we /jave 

ll^jWnll <4|t;| 1 (ri+l) 1/2 ||^|| 

Proof. By the choice 2 J > n we do not need the projection operators Pj, in other words 
V n C 7"0 = PjT . For simplicity, denote 

Ae = v~gj{x) 1/2 a x 

We have 

F J ( w ) = ^A; + ^-(^) 1/2 ^ aX 

X 

and 

||FjWn|| < ^||K+^ n || + (jY /2 J2\ d M< a ^n\\ 

x x 

The second term is bounded by 2 1/2 J -1/2 n|i;|i||^> TI ,||; f° r the first term we use 
||(4; + A x )^ n f < 2\\(A:.A x ) 1/2 ^ n \\ 2 + 2\\{A x A* x y^ n f 
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On T> n , we have, using 2 J > n as well, 



A X A X 



\v x \ 2 a* x gj(x)a x = \v x \ 2 n ^( 1 ~^j) <\v x \ 2 {n+l)t 



A x Al = \v x \ 2 gj(x) 1/2 a x a* x gj(x) 1/2 =\v x \ 2 (n x + l)(l-^) < |^| 2 (n + l)l 



and, using once more 2 J > n 
\\Fj(v)i/) n \\ < 2\v 



1 Wn 



(n + l) 1/2 + 



?? 



(2J) 



1/2 



< 4\v U(n + 1 



|V2| 



□ 



Denote with s-lim the strong resolvent operator limit for operators acting on JF. The 
following is our general result about the convergence of the fluctuation observables. 

Theorem 4.4. For all k G N and i>x> • • • > u fc £ ^;s(£)> u ' e ^ aue 

s-hmn^(^)=n F ^) 

where the limit is taken over any sequence of J G |No tending to oo. More precisely, for 
all n G N, ?/> n G V n , and J G 4N suc/i i/iat 2J > n + k, we have 



k k 
3=1 3=1 



< 



(2jy 



/2 



k 

n 



«i i 



(n + fc)hV2 



n! 



^ 2 fc +> + £;-^ 2 ||^| 



(9) 



It is convenient to single out the k = 1 CclSG clS cl separate result, and prove this first, 

which we do in the following lemma. The proof of the Theorem 14.41 is then given after 
the proof of Lemma 14.51 

Lemma 4.5. For v G we have 



s-lim Fj(v) = F(v) 



(10) 



where the limit is taken over any sequence of J G |No tending to oo. More precisely, for 
all n G N, ip n G T> n , and J G |No such that 2 J > n + 1, we have 

n 



[F{v)-Fj{v)]^ n \\<Mi^jy^\\^ 



(11) 



Proof. First note that (jl(Jj) follows directly from (Jllj) . see, e.g., Theorem VI 1 1.25 of [21], so 
it is sufficient to prove This can be done as in the proof of Lemma [4.31 This time, 

A x = v x (l- gj(x) 1/2 )a x 

We have 

|| [F(v) - Fj(v)]^ n \\ < Y,\\( A t + A Mn\\ + (§) 1/2 E I^HhX^H 



CENTRAL LIMIT THEOREMS FOR THE ASYMPTOTICS OF QUANTUM SPINS 

The second term is bounded by 2 1/2 J _1/2 n|u|i||^ n ||; for the first term we use again 

||(^ + A^n|r<2||(A:A,) 1/ Vn|r + 2||(A x A:) 1/ Vr l |r 
These two terms can be bounded using 1 — gj(x) 1/2 < n x (2J)~ 1 and Lemma f4. 21 

\\{A* x A x ) 1/2 i) n \\ 2 = (V> n ,A*A^n> = \v x \ 2 (a x ^ n ,(l - g.j(x) 1/2 ) 2 a x ij n ) 



13 



. w , q y f t \ 2 I I iiO I I 9 / ^ \ 2 

< 1^1 [—) \\a x ip n \\ < \v x \ {—) v\ 



A x A x ) 1/2/ ip n \\ = (i> n ,A x A* x i; n ) = \d x \ 2 (^ n ,(l - gj (x) 1/2 )a x al(l - gj(x) 1/2 ^ n ) 
<\Vz\ 2 (^) 2 (n + l)\\M 2 



Putting everything together, and using 2 J > n + 1, we find 



\\[F{v) - Fj{v)]^ n \\ <2\v\ 



n(n + 1) 1/2 ^ n 



23 



(2J) 



1/2 



.11 <4M 



7? 



(2J) 



1/2 



□ 



Now follows the remainder of the proof of Theorem 14.41 



Proof of Theorem As in the proof of Lemma 14.51 it is again sufficient to prove eq. 
(jHJ). We write the difference of products as a telescopic sum: 

k k 

Y[F(v s ) -HFj( Vj ) = [Ffa) - Fj( Vl )]F(v 2 ) . ..F(v k ) 

3=1 3=1 

+ Fj( Vl ) [F(v 2 ) - Fj(v 2 )] F(v 3 ) . . . F(v k ) 
+ ... 

+ Fj(vi) . . . Fj{v k „ x) [F(v k ) - Fj(v k )] 

We estimate each term separately, first using repeatedly Lemma [4. 3[ then using Lemma 
14. 5\ and finally using repeatedly Lemma f4. 21 

Fj(«i) . . . F^-i) [F(vi) - Fj(vi)]F(v i+1 ) . . . F(v k )i/> n 

i(n + k — j + l) 1 ^ || [Ffr) - Fj(vA] F(v i+1 ) . . . F(4)Vv, 

(n + k — i) 



i-l 



< 



< 



[ri 4 M 

i-l 

YlAlv^n + k-j + l) 1 ^ 4| Vi | r 



i-l 



(2J) 



1/2 



< [n 4 i^i i ( n+A; -^ +i ) i i 4 i^i i 

3=1 



(n + k — i) 
(2J)V2 



. . . F(v k )ip ri 



J] 2|^| 2 (n + A;-i + l) 



1/2 
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< 



[n 



Vj\l 



ft 

[n 

i=i 



4*2 



(2J) 1 ' 2 
(n + fe)! (n + fc 



li[(n + fc-i + lH||V» 



n! 



2J 



V2, 



Vv, 



□ 



Now we return to the main subject of this paper, namely proving central limit theorems 
for the coherent state us. 

We get convergence of all moments as a trivial application of the Cauchy-Schwarz 
inequality and Theorem 14 A\ more specifically eq. © for the case n = 0. 

Corollary 4.6 (Moments). For all k G N and t>i, . . . G ^(-C); we have 



lim a; 

J^oo 



A V 3 



CO 



3=1 



where the limit is taken over any sequence of J G ^N . More precisely, for all J G 4N 



— a; 



< 



2' 

(k\) 1/2 



(2J) 



1/2 



k 

n 



1/2 



It is clear that Theorem 14.41 contains much more information about the convergence 
of the fluctuation operators than the convergence of all moments that is derived from 
it, or the convergence of characteristic functions that is given in Theorem 12 .HI In the 
usual setting of quantum central limit theorems [111101 a space on which all fluctuation 
operators as well as the limiting boson field operator act simultaneaously, does not exist 
and therefore the question of strong (or any other) operator convergence does not make 
sense in these situations. It is an interesting question, however, to ask whether the 
existing central limit theorems in the usual setting for sums of random variables can be 
strengthened. Therefore, we end this section with several reformulations of Theorem 14.41 
to obtain statements that do make sense although they have not been proved in general. 
In the Section IB~2l we will give a first example of such a stronger convergence result. 

Corollary 4.7 (Spectral measure). For any k G N and v\,...,Vk G £j(£), the spectral 
measure in the coherent state u of Fj(v\) . . .Fj(v)~) converges to the spectral measure in 
the Fock state uj of F{v\) . . . F(v k ) , where convergence is as functionals on the 
bounded continuous functions on R. 

Proof. It follows from Theorem 14.41 and Theorem VI 1.20 of [21] that for any / G 
and any if) G T 

k k 



jiv(n^)>=/(n^ 



■0 
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□ 

The following reformulation is probably the most interesting. It generalizes Theorem 
12. HI by proving convergence of the characteristic function of arbitrary polynomials of fluc- 
tuation operators instead of only allowing linear combinations of them. 

If A\ , . . . , Ak are selfadjoint elements of a general C* or von Neumann algebra 21, we de- 
note with C(Ai, . . . , Ak) the ring of noncommutative polynomials in Ai, ... , A k . This ring 
admits a *-involution which conjugates each coefficient and reverses the order of multipli- 
cation in each term. A polynomial p G C(A\, . . . , Ak) is called a selfadjoint polynomial if 
it is invariant under this involution. For example, the anticommutator A\A 2 + A 2 Ai is a 
selfadjoint polynomial, but the monomial A\A 2 is not. 

Corollary 4.8 (Characteristic function of polynomials). For any k G N, Vi,...,Vk G 
and p G C(Fj(vi), . . . , Fj(v k )) a selfadjoint polynomial in k variables, we have 

lim uj(y p[Fj( - vl) '---' FjM] ^ = u (ypl F (vi)>---,F(v k )]^ 

where the limit is taken over any sequence J G |No tending to oo. 
Proof. It follows from Theorem 14.41 that 

s-limp^i), . . . , Fj{v k )\ = plFfa), . . . , F(v k )\ 

and by Trotter's theorem (Theorem VI 1 1.21 of |21j ) also 

s-lime ip[Fj{vi) '-' Fj(vk)] = e ^(«i),-,-F(fife)] 

J— >oo 

The result follows from the Cauchy-Schwarz inequality in T . □ 

Remark 4.9. In the previous corollary we could further generalize and consider I selfad- 
joint polynomials pj in kj variables. It follows by the same argument that 

Jim - (n^" 7 :/i< : ' /Vr v) = &(f[e ipjlF{ ^ 

3=1 3=1 

In addition all products of polynomials and characteristic functions of polynomials (e.g., 
Pi x e tp2 ) can be considered as well. 

Remark 4.10. We have always taken v G ^(-C). This is convenient because then the 
finite J fluctuation operators Fj(v) are bounded operators with 

1 1 || < 2 3/ >| x J 1/2 

We can ask about more general examples where Fj(v) is only densely defined, and one 
such generalization is the following. Consider in expression (jHJ) the last term 

Vj) l^ V * a x ax 
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On V n this operator is bounded by 2 1/2 J 1/8 n sup xe£ \v^\- Hence all of our results remain 
unchanged, except for notationally more complicated error bounds, if we consider v E 
£3(2) instead of v E ^(-C), where £3(2) is defined as the set of Revalued sequences on £ 
such that v = {(^,^)} x6£ e ^(£) and £ 3 = {{g} xe c e £f{T). 

Remark 4.11. In all of the previous results, we can replace the state Qj = Q by a 
perturbed state flp, as long as the perturbed state is still in the set D of analytic vectors 
for F(v). 

5. Applications 

5.1. Bosonization for quantum spin Hamiltonians. The main application of Theo- 
rem 14.41 is for the large-spin asymptotics of quantum spin Hamiltonians and their corre- 
sponding dynamics. Consider for instance an interaction between spins in a finite volume 
A of the type 

3 3 
H ^ = "EE h ij (x,y)S i x Si - J£;X>(4S£ (12) 

x,y£A i,j=l xgA i=l 

where the interaction functions and g>j satisfy all necessary conditions for selfadjointness 
of -£/j,a, and in addition are bounded and of short enough range, i.e., 

su.p(y2\h i:j (x,y)\j < 00 
sup \gi(x)\ < 00 

More general Hamiltonians can easily be treated along the same lines, but this example 
already includes the various interesting Heisenberg-type models. The minus sign in front 
of the Hamiltonian is there for convenience, we do not make any explicit assumptions on 
the sign of the hij or g>j. 

Bosonization refers to the idea that the low-energy excitations of a quantum spin system 
can be effectively described by a boson approximation to ()12j) . describing non-interacting 
bosons hopping on the lattice £. To make this into a mathematically precise statement, 
the Hamiltonian has to be scaled by J~ l and a large-spin limit has to be taken. Again there 
exist various ways of taking this limit. One way consists of proving that the free energy at 
low temperatures converges to the free energy of a boson model, see, e.g. |4" Ho1l22| . and this 
has been used mainly to study the isotropic Heisenberg ferro magnet. Another approach 
was introduced in ^H] i n the study of 1-dimensional anisotropic Heisenberg ferromagnet, 
and consists of proving strong operator convergence of the Hamiltonian like we did for 
the fluctuation observables in Theorem 14.41 From there, convergence of eigenvalues and 
eigenvectors is derived. 

In both approaches it is necessary that there is a unique ground state to be used 
as a reference state, and this ground state should be a coherent product state. Both 
the isotropic and anisotropic Heisenberg model have a rotational symmetry, and hence 
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uniqueness of the ground state can only be achieved by adding an external field which 
breaks this symmetry. For the isotropic model, there is full rotational symmetry, and 
the final results do not depend on which ground state is selected. For the anisotropic 
model, there exist non translation invariant ground states describing interfaces, and the 
approximating boson model is different for different ground states, see [IB] for details. 

In this section we want to demonstrate another use of Theorem 14 A\ namely in studying 
the dynamics of (ground) states under perturbations. The main assumption is that there 
exists a ground state uj of (fT2|) which is a product state of coherent states, like we used 
throughout this paper. For the present application it is not necessary that it is the unique 
ground state. 

The first step is to write the Hamiltonian (|12|) in the new natural basis for the state oj: 

3 3 
x,y£Ai,j=l x£A i=l 

where the hij and g~i are easily obtained using the basis transformation on page|Hl Instead 
of taking i,j = 1, 2, 3 we can also take i,j = — , +, 3. The coherent state in the new basis 
is just the 'all up' state, and this is an eigenstate of Hj\ if and only if the only non-zero 
coefficients are as follows: 

Hj,a = ~ E [h-+{ x ,y)Sx~Sy + ~h+-(x,y)S+Sy + h 33 (x,y)S x Sy - J^2g 3 (x)S x 

x,y£A xGA 

- E [^(^ y)s x Sy + hi{x, y)SxSl + ^sO, y)s 2 x sl + h 32 (x, y)S*S*\ 

x,yeA 

For the terms on the first line, no additional conditions except for those of selfadjointness 
and summability are needed, the terms on the second line however can only be allowed if 
for all x G A 

E ^ 3 ( x ' y) = E y"> = for i = X ' 2 

yeA yeA 

Suppose the reference state u is a ground state of Hj\. Then, its energy is given by 

x,y&A xdA 

Hence we can renormalize the Hamiltonian such that the ground state energy is 0, and 
obtain 

Hj, A = -£ [h-+(x,y)§-§i + ~h + -(xiy)~St~Sy + h3 3 (x,y){S 3 x S 3 y - J 2 ) 

x,y£A 



-jJ2h(x)(S 3 x -J) 



zga 
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13 (x,y)Sl(S 3 y - J) + h 31 (x,y)(S 3 x - J)^ 1 



+ h 23 (x,y)S 2 x (S 3 y - J) + h 32 (x,y)(S 3 x - j)^ 

With the Hamiltonian correctly renormalized, we can also write down the GNS Hamil- 
tonian Hj acting on the GNS Hilbert space Tij, it is the same expression with all sums 
extended from A to the entire lattice £. 

It is also clear that after scaling with J -1 , this Hamiltonian is of the form 

^=-EpSW)'^.^))4 E ~h 3 3(x,y){J 2 -SlS 3 y )+J2h(x){J-S 3 x ) 



cc,ye£ 



x£A 



(2) 

where Pxy are selfadjoint polynomials in 2 of the 3 listed variables, and /* are the rotated 
basis vectors defined on page |H1 

Using Theorem 14.41 we get immediately that 



x,y££ x,ydA 

The remaining terms can be evaluated by writing them on Fock space: 

E 9^ x ) { J - ^) = E 93(x)n x 

x£A 



x&A 



and 



j E M*,y)(J 2 -^ 3 ) = E m*,j/) 



J - j(J - n x )(J - n y ) 

x,y££ 

^ ^33 (x, y){n x + n y - jn x n y ) 

x,y€£ 

^ E(x)n x — h^ix^) 7 ^ 1 - 

xe£ a;,J/e£ 



where 



= E Pssfo + ^ 33 (^> 

Using Ylx n x = on it is seen that the second term above converges strongly to 0. 
Hence we conclude: 

Sjlim ±Hj = H := ^2[E(x) + g 3 (x)]a* x a x - ^ [h-+(x, y)a* x a y + h + _(x,y)a x a* y ] (13) 

i££ x,y&A 

and this is a rigorous way to express bosonization in a quantum spin model which is 
of the type (fEZj) and has a coherent ground state. As the coefficients g 3 (x), h_ + (x,y), 
etc., depend on the rotation of the basis, it is clear from this derivation that the boson 
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Hamiltonian in general will be different for different ground states. Notice also that it is 
the second quantization of an operator H on £ 2 (£) defined by: 

(Hv) x = [E(x) + g 3 (x)]v x -J2[h-+( x ,y) + h+-(y, x )] v y ( 14 ) 

In [TH], eq. (|T3j) is used to show that the eigenvalues and eigenvectors of the boson 
Hamiltonian give the large-spin asymptotics of the eigenvalues and eigenvectors of the 
spin Hamiltonian (for the special case of the anisotropic ferromagnetic Heisenberg chain). 
This is a very practical result because the spectrum of the spin Hamiltonian is generally 
very hard, if not impossible, to compute for large systems, while the spectrum of the 
boson Hamiltonian, which is quadratic, can be computed very easily, also for large systems 
(certainly numerically). 

Theorem 14.41 can also be used to obtain a rigorous estimate for another kind of compu- 
tation, namely for the dynamics of states, typically eigenstates of Hj, under perturbations 
of the Hamiltonian. Such a problem is typically studied in a weak coupling limit, where 
the strength of the perturbation vanishes as the microscopic time tends to infinity, see 
[TT?] for a recent example studying the dynamics of interfaces. 

Clearly, if we perturb the spin Hamiltonian by any selfadjoint polynomial in the fluc- 
tuation operators, we can study the influence of this perturbation using the bosonization 
approximation, i.e., taking the large-spin limit. More precisely, let p be a selfadjoint 
polynomial in k variables and let v\, . . . , V}~ G P\{£)- Denote 

Pj = p[Fj( Vl ), Fj(v k )] P = p[F(vi), F{v k )] 

Then 

s-lim ±Hj + Pj = H + P 

J— >oo J 

and again by Trotter's theorem (Theorem VIII. 21 of |21j). for all ip G J 7 , 

lim e^^+^V = e^ +P V 

J^oo 

While the spin expression on the l.h.s. can again be very difficult to compute, the boson 
expression on the r.h.s. can be solved for relevant cases. 

Consider for instance the case where the Hamiltonian is perturbed by a fluctuation of 
the spin in a certain direction, i.e., Pj = Fj(v) for some v G £\(£)- Such a perturbation 
was also considered in in the study of spatial fluctuations of quantum spin systems 
and their relation with linear response theory. Using a Dyson expansion and the CCR 
algebraic structure we find 

e it[H+F(v)] e -itH = e -±fids<r(v,v') e iF(&) 

where 

v*= f dse lsH v 
Jo 
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and H is given by ()14j) . Hence the boson approximation to the time evolution of the 
ground state under the perturbed evolution is given by 

Using techniques such as in the previous section and such as in Lemma rA.2l below. explicit 
estimates which vanish for large values of J can be obtained to compare this state with 
the true evolved state 

e it(iH J+ p. 7 ) fij 

5.2. Strong central limit theorem for the large N limit of iV spin-| particles. In 

this section, we show that Corollary 14 . 81 gives a fully noncommutative example of Theorem 
1 of 14 j. In this example, we consider a limit where the number of spin variables, N, 
tends to infinity, while the magnitude of each spin remains constant equal to 1/2. We 
change the notation accordingly. 

A quantum probability space consists of a C* or von Neumann algebra 21, and a state 
lu on 21. Then (%l® N ,U]y = u® N ) denotes N independent copies of (21, u), and if A G 21, 
then 

Aj = l<g)---<g)l<g)^4<g)l<g)---<g)l 

with A in the i th position, denotes the i th copy of A. 

The fluctuation operator of A G 21 in this context is defined as 

1 N 

F N (A) = - / =J2^- UJ ( A )] 

* i=l 

and we can ask the same convergence questions as we did before. In fact, this simple 
setup was the first in which the connection between quantum fluctuation operators and 
representations of the CCR was worked out in detail [7j. In that paper, the result analogous 
to Theorem 12.31 was first proved. 

The question whether stronger central limit theorems holds in this setup was first raised 
in and in the following special case was proved: 

Theorem (Kuperberg [T4"]). Let (21, r) be a quantum probability space with r a tracial 
state. For allk G N, A\, . . . , A^ selfadjoint elements o/2l, andp G C(Fjv(Al), . . . , F^(Ak)) 
a selfadjoint polynomial in k variables, we have 

lim rw (V p[Fiv(Al) '-' Fjv(j4fc)] ) = E ^ e i P [x(A 1 ),...,x(A k )]^ 

where X(Ai), . . . ,X(Ak) are classical Gaussian random variables with covariance matrix 

E(X(A i )X(A j )) = r(A t A 3 ) 

The fact that a tracial state leads to classical fluctuations in the limit is to be expected, 
it follows directly from the cyclicity of the trace. As far as we know, no generalization of 
this theorem to non-tracial states exists. 
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Consider the simple case in which 21 = M 2 (C), the complex (2 x 2)-matrices. 21 carries 
the irreducible spin-| representation of 57/(2) given by | times the Pauli matrices 

*' = (° D ^ = (° i\ o'-( l 



i oy v / v° _1 , 

The standard basis of C 2 is diagonal for cr 3 , we denote |1) = (J) and | - 1) = g). For the 
state we take u> = ujg t!fi , the expectation in the spin-| coherent state 

|(^)> = e ^^-- +e ^)|l) 

where o~ ± = a 1 ± io 2 . 

We are interested in the selfadjoint operators 

v a = uV 1 + v 2 a 2 + v 3 a 3 , v G R 3 

and their fluctuation operators Fn(v) := F^(v ■ a). 

But a collection of N independent spin-| degrees of freedom can also be regarded as 
one spin-y degree of freedom, i.e., define 



N 

J 

i=l 



then these operators satisfy the SU(2) commutators as well, and the eigenvalues of each 
Si obviously are — — + 1, . . . , — . Moreover we have in the coherent state to: 



N 



i=l 

s(§M) + i-r-r C0S (§M) 



{ 



N 




where u = (9,tp). This is the correct generating function for the spin— coherent state, 
cfr. eq. ((21). Hence the present setup is exactly equivalent with the setup in the pre- 
vious sections if we take the lattice £ consisting of a single point, and check that the 
normalization of the fluctuation operators is consistent with eq. (j2J): 

2 N 

[v ■ S N - u N (v ■ S N )] = —= V] [\v ■ Ui - \u(v ■ (Ti)] = F N (v) 

It follows that Corollary 14. 81 remains valid in the present setup, and we get immediately 
the follwoing theorem. 

Theorem 5.1. Consider the quantum probability space (IK^C), to) where uj is the coherent 
state defined by the unit vector u = (6,(p) G R 3 . For all k G IN, V\,...,Vk G R 3 , and 
p G C{Fn(vi), . . . , F]y(vk)) a selfadjoint polynomial in k variables, we have 

lim UN (jPlFN(vi),.,F N (v k )]\ = ~( e i P [F(v 1 ),...,F(v k )] 
N^oo V / V 
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where u is the Fock state on CCR(IR 2 , a), a is the standard symplectic form, and F(v) are 
the boson field operators for tie R 2 which is obtained from v e R 3 by projecting onto the 
tangent plane to the unit sphere at u. 

Note that the limiting object in the above CLT, the Fock state on CCR([R 2 ,o~), is es- 
sentially the quantum harmonic oscillator. It can be regarded as a non-commuting pair 
( "position" and "momentum" ) of Gaussian random variables. 



Appendix A. Proofs for Section [21 



Proof of Proposition WHl The characteristic function of Fj(v) is, using eq. (J2J) and some 
straightforward manipulation, 



UJ 



( e *HM«)) = Y[e' it ^ uj{Vx ' Sx) uj(e it ^ Vx ' Sx ) 



From the formula 



n 



-it{2J) 1 / 2 v x -u 1 



COS 



Vx ■ U 3 
\Vx\ 



t\Vx\ s 

(2J) 1 /2- 
it\v. 



+ i- 



.V x ■ Un 



sin 



t\v x \ 



cxp 



CO 



1 A v x - u x 



(e^W) = exp{^ ^ T ( Fj(v), Fj( v) 
k=l ' 



(2J) 1 / 2 
it\v r 



\v x \ *(2J)v*' 



2.1 



(2jy 



/2 



1 + 



\ V x\ 
Vn- ■ U 



\v T \ 



2.1 



it follows that the fc-point truncated correlation function is given by 



t=0 



t=0 



where f x (t) is the part between {. . . } above. Clearly f x (0) = 1, f' x (0) = 0, and hence 

(w,)'(o) = (mo(o) = o 

Ob/.Ho) = (m - C£)7£)(o) = f'M 



(ln/,) (fc) (0) = /i fc) (0) 
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We find for k>2: 
lo t ( Fj{v),...,Fj(v)) 
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k 



+ 2 (1 



v x • u. 



(2J) 



i\v 3 



1/2 



(2J) 



1/2 



1 - 



1 + 



\V X \ 



\v 



■I' I 



l k 



- - ^){ (i - + «- 1 > t ( 1 + ^r)1 

For /c = 2 this reduces to eq. i.e., u;t(-Fj(w), Fj(v)) = \vfy, while for > 2, we get an 
upper bound 



such that 



OO .J, 



fc=3 



fc! 



fc=3 



Jl/2 



We have 



o;( e ^(«)) -Cj(W(v))\ 

1 °° 1 

< |exp{--|j?| 2 + ^2^ut{ Fj(v),...,Fj(v)) ] -expj-- 



k=3 



k 



oo 

^ l ex p(E ^t( Fj(v), Fj(v)) | 

fe=3 ' ,„ 



For z e C 



and hence 



e z - 1| = I / dse sz z| < e lzl \z\ 



uj(e iFj( - v) ) -u(W(v)) 



2 1/ >| 2 + 2 1/2 exp(2 1/ >| 2 ) 



□ 



where we used J > ~ in the exponent on the r.h.s. 

Proof of Proposition \2.1\ The proof follows from the previous proof. Consider the centred 
observable 

1 x 1 

M j(v) = - ^2(v x ■ S x - J{V X ■ U x )) = ( 2 jy /2 F j( V ) 
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u t { Mj(v),...,Mj(v)) = ——cu t ( Fj(v),.. : ,Fj(v)) 



and from eq. (fT5j) . for k > 2: 



u t {Mj{v),...,Mj{v)) 



< 



Mi 



-i 



Continuing as in the proof of Proposition 12.21 we find 



< 



(0 



( e "M«)) _ 1 



< 7 exp 



Mi 



+ 2exp(|u|i) 



The first step in the proof of Theorem 12.31 is the following Lemma. 
Lemma A.l (A BCH-formula). Let Vi,v 2 G t/ien, /or J G |N , 



giFj(Di) e jFj(D 2 ) _ giF J (t;i+ i i2)g-|[Fj(Di),Fj(j;2)] 



< 



1 



3J 1 / 2 



lfa.|2| v 2 | 2 (|f a.| 2 + |u 2 | 2 ) 



where \\ ■ \\ is the operator norm in B(7ij). 

This Lemma is a special case of 
Lemma A. 2. Let 21 be a C* -algebra, A, G 21 selfadjoint, and C G 21 arbitrary, then 



e^\ L7]||<|p,L7]|| 



iAJB _ i(A+B) — | [A,B] 



e e 



<^[\\[A,[A,B))\\ + \\[B,[B,A))\ 



Proof. Eq. (jTB| is well-known and proved as follows: 

gfaA^g-faA = c + i f dt e i{s - t)A [A, Cle-* 
Jo 

and hence 



(s-t)A 



JsA 



,C\=i / dt e^'- t)A [A, C\e { 



itA 



Therefore 



Next we prove (fTTj) . 



[e iA ,L7]|| < / dt \\e^ l - t)A [A,C]e itA \\ < \\[A,C]\\ 



e *V*(l _ g-iBg-iA ei (A +J B)g-|[A, B ]) 

< 1 _ e -iB e -U -I[A,B] 



e iA e iB - e i(A+B) e -±{A,B] 



(1 
(1 



(1 
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Let 

F{t) = 1 - e -iB e -UA e i(tA+B) e -HA,B] 
We compute F'(t): we have 

diitA+B) = • f 1 ds i.(tA+B) ^ i(l-.)(tA+B) (lg) 

dt J K J 

and hence 

F'(t) = te- w e- UA Ae^ tA+B ^e-^ - i f ds e~ iB e -^ e is(tA+B) Ae i { i- s){ tA+B) e -|L4,B] 

Jo 

+ ^ e -iB e -itA e i(tA+B)^ B ^i[A,B] 

Compute, using (JTSJ), 

lAe i(tA+B) _ • / dfj e is(tA + B) Ae i(l-s)(tA+B) 



and also 



= i f 1 ds (Ae is{tA+B) e i{1 ~ s){tA+B) - e is{tA+B) Ae i{l ~ s){tA+B) ) 
Jo 

= i f ds [ A ^ e ^tA + B) ]e i { l-s)(tA + B) 

Jo 

= - f ds f du e Ks-u)(tA+B)^ B y ( i- s +u)(tA+B) 

Jo Jo 

l e i(tA+B) ^ S ] = f ds f du e i( s -u)(tA + B) e i(l-s + u)(tA+B) ^ ^ 

2 Jo Jo 



Hence 

-1 r s 



We have 



F'(t)= [ ds rrf M e-^e- i4 V (s -" )(M+B) [e i(1 - s+ " )(M+B \[A5]]e~5^ 
Jo Jo 

F(l) = F(0) + [ dt F\t) 
Jo 



and hence, using (fTb]) . 

-1 r i 



\\F(1)\\< [ dt f ds f du [e i{1 - s+u)(tA+B \[A,B]] 
Jo Jo Jo 

/*1 pi PS 

< dt ds du\l-s + u\\\[tA + B,[A,B]]\\ 
Jo Jo Jo 

< J ds J du\l - s + u\ (\\[A, [A, B}}\\ + \\[B, [A, B]]\ 
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e iA e iB _ e i(A+B) e -\[A,B] 



<-(\\[A^B]]\\ + \\[B,[B,A))\ 



□ 



Proof of Lemma \A.l\ First note that by the SU(2) commutation relations 

1 x 1 . 

[Fj(v), Fj(w)] = -f^lvx- S x , w x ■ S x ] = — y^ y ( v x x w x) ' Sx 



and hence by the previous lemma 



e iFj(v) e iFj(w) _ e iFj(v+w) e -\[Fj{v),Fj{w) 



< -^j^ ^|| (v x x (v x X w x ) + w x x (w x X v x )) ■ S a 



- x ( Vx x Wx ^ +w * x ( w x x 



< t^7 2 Yl \ v *\\ w *\(\ v *\ + KD ^ ^T^l^bkbdf | 2 + k| 2 ) 



3J 1/2 



□ 



Proof of Theorem \2.tA It remains to prove the induction step. Assume (jHJ) holds for some 
n G N , and choose v±, . . . ,v n+ i G ^1(£). For convenience, denote 



By the BCH-formula 



u)(Wj\e iFj ^e iFj{ - Vn+1 ^ — e iFj ( Vn+Vn+1 ^e~^ Fj ( Vn ^ Fj( - Vn+1 ^~\ J 



< 



3J 1 / 2 

Using the Cauchy-Schwarz inequality, and Proposition 12.11 



±[F.j(v n ),Fj(v n+1 )]\ e -±a(v n ,v n+1 ) 



uj\e 2 



\[Fj(v n ),Fj(v n +i)]\ e la{v n ,v n+1 ) 



< — exp 

•J 



\v n \2\V n +l\2 + 2 exp ( |f n 1 2 K+l 1 2) 
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Using these two inequalities, and the CCR-algebraic structure, we find 
(w je iFj(v n ) e iFj(v n+1 )^ _ fi^ W ( Vl ) . . . W(vSj 



< 
+ 
+ 



W je i F j(vn) e iFj(v n+1 )^j _ u ^y je iFj(vn+v n+1 ) e -±[Fj(v n ),Fj(vn+i)]^ j 

uj (Wje iFj ( Vn+Vn+1 ^ e~^ Fj ( Vn ^' Fj( > Vn+1 ^] j - u(w J e iFj( - Vn+Vn+1 ^ei' 7 

tJ[Wj ^Fj(v n +v n+1 )\^(v n ,v n+1 ) _ . . . W(v n -!)W(v n + V n+1 ) ) e ~^,v n+1 ) 



r(vn,v n +l) 



u W{v x ) . . . W{v n ^)W{v n + v n+1 ) e 



'i<r(v n ,v n+ i) _ ~ 



Q{W(v 1 )...W(v n ) 



< 



+ 



Wje iFj(v n+ v n+1 )^ _ Q [ W ( v j W(v n ^)W(v n + u n+1 )) 



UJ 

i r l 



, i . JI -'n|2|Un+l|2(|Un|2 + |Un+l|2) + V^exp - | W n | 2 | V n +1 b + exp ( | V n \ 2 \ V n+1 1 2 ) 



In the second term, the part between {. . . } is the function a(v n ,v n+ i) as defined in the 
statement of the theorem. Now use the induction hypothesis on the first term, and the 
desired result follows. □ 
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